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RESEARCH MEMORANDUM 


THEORETICAL AND EXPERIMENTAL INVESTIGATION OF 
ARBITRARY ASPECT RATIO, SUPERCAVITATING 
HYDROFOILS OPERATING NEAR THE 
FREE WATER SURFACE 
By Virgil E. Johnson, Jr. 

SUMMARY 


Theoretical expressions for the lift, drag, and center of pressure 
of supercavitating hydrofoils of arbitrary section, aspect ratio, and 
depth of submersion are developed for the zero-cavitation-number case. 

In order to verify and supplement the theoretical investigation, 
an experimental investigation was made in Langley tank no. 2 of two 
aspect-ratio-1 hydrofoils, one with a flat and one with a cambered lower 
surface. Zero cavitation number was obtained in the tank by operating 
the hydrofoils near the water surface so that their upper surfaces were 
completely ventilated. Data were also obtained at cavitation numbers 
greater than zero for flows with vapor-filled cavities and with air arti- 
ficially introduced to simulate the same flows at lower speeds. The tests 
covered a range of water speeds up to 80 feet per second, angles of 
attack up to 20°, and depths of submersion from 0- to 85-percent chord. 

The theoretical force and moment coefficients agreed with the data 
obtained in the zero-cavitation-number tests within an accuracy of about 
5 percent. The ventilated force and moment coefficients of both models 
were independent of speed. As predicted by the linear theory of Tulin 
and Burkart, the cambered hydrofoil lift-drag ratio was superior to that 
for the flat plate at the same lift coefficient. Similar to the influ- 
ence of camber on a fully wetted airfoil, the influence of camber in 
ventilated flow was to effectively increase the angle of attack without 
appreciably changing the lift-curve slope. Comparison of the present 
results with the results of other investigations revealed that a sharp 
leading edge is necessary for good efficiency in supercavitating flow. 
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The nonzero-cavitation-number tests confirmed the finding of others 
that the lift of a hydrofoil operating in a finite cavity depends only 
on the cavitation number and is independent of whether the gas in the 
cavity is air or vapor. 


INTRODUCTION 


The desirability of using auxiliary lifting surfaces such as hydro- 
foils or hydroskis for reducing seaplane-hull loads and improving rough- 
water performance is well established. Although the higher aspect-ratio 
submerged hydrofoils with tjieir accompanying high efficiencies could be 
the most desirable lifting surfaces to be used, only the low-aspect-ratio 
planing hydroski has so far been successfully applied as landing gear to 
modern high-speed aircraft. The reason for this is that the conventional 
hydrofoil presents problems not experienced by a hydroski. 

As a hydrofoil nears the free water surface (during a take-off run) 
the low pressure side of the hydrofoil becomes ventilated from the atmos- 
phere. This phenomenon results in a severe and usually abrupt loss in 
lift and reduction in the lift-drag ratio. For conventional airfoil 
sections, the loss in lift may exceed 75 percent. The speed at the 
inception of ventilation depends on the angle of attack and depth of 
submersion; but, except for very small angles of attack and relatively 
low take-off speeds, the speed at inception is usually well below the 
take-off speed of the aircraft. 

Even if the ventils Mon problem is overcome by using small angles 
of attack and incorporating fences or other devices for suppressing 
ventilation, the onset of cavitation presents a second deterrent to the 
use of conventional hydrofoils at high speeds. The loss in lift accom- 
panying cavitation of conventional airfoil sections is not abrupt, but 
the ultimate reduction in lift and lift-drag ratio is comparable to that 
of ventilated flow. Even thin-airfoil sections of small design-lift 
coefficient enter this cavitating regime of poor lift-drag ratios at 
speeds in excess of about 80 knots. 

Since the take-off speed of modern supersonic aircraft is in the 
range of 150 to 200 knots, lifting surfaces with cavitating or venti- 
lating characteristics superior to those of conventional airfoil sec- 
tions are desirable. Fortunately, theoretical work by Tulin and Burkart 
(ref. 1) has shown that superior configurations do exist and they have 
selected a cambered configuration for operation in cavitating or venti- 
lated flow which, in an inviscid fluid, has two-dimensional lift-drag 
ratios at its design angle of attack and zero cavitation number, about 
six times that of a flat plate. If such a cambered hydrofoil can be 
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induced to ventilate at very low speeds, while the aircraft hull still 
supports most of the load, a stable and efficient take-off run may he 
possible. 

The present report is concerned with an experimental and theoretical 
investigation of one of the Tulin-Burkart cambered configurations and a 
comparison of its characteristics with those of a flat plate. The pur- 
pose of the investigation was (l) to find out if the superiority of the 
cambered configurations obtained from two-dimensional, infinite -depth 
theory would be maintained at finite depth and aspect ratio, and (2) to 
develop expressions for the prediction of the force and moment character- 
istics of a ventilated hydrofoil of arbitrary section profile which 
would include these three-dimensional effects. 



SYMBOLS 

A 

aspect ratio 

Anj-^O 

coefficients of sine-series expansion of vorticity distri- 
bution on equivalent airfoil section at infinite depth, 
that is, 

ft(x) = 2V (aq cot ^ + Aj_ sin 0 + A2 sin 20 . . . A n sin n0^ 
where 

A n = — /^ cos n0 d0 

* J 0 to 

A/~\ = — — /”* — ^ d 0 + (X = (X + A/~\ 1 
^ n J 0 cbc ^ 

An,h>Ao,h 

coefficients of sine-series expansion of vorticity distribu- 
tion on hydrofoil section 

a 

distance from equivalent airfoil leading edge to center of 
pressure in chords 

Bn^O 

coefficients of cosine series defining location of image 
vortex in airfoil plane 

b 

parameter defining location of spray at infinity in t plane 
(see ref. 2) 
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'D 


C D,1 


■'DjC 


'L,l 


'L,c 


- / m 


total-drag coefficient, — 

qS 

profile-drag coefficient (total drag less crossflow and 
. Dn 

skin friction), — 

qs 

crossflow-drag coefficient (normal to the section) , 

Df 

skin-friction drag coefficient, — 

' qS 


total-lift coefficient, 


qS 


total- lift coefficient of equivalent airfoil section, — 

is 

Li 

lift coefficient exclusive of crossflow, — 

qS 


crossflow lift coefficient, — £■ 

qS 

pitching-moment coefficient (about the leading edge), — — 

qSc 


-m 


pitching-moment coefficient of equivalent airfoil section 
(about the leading edge), 

qSc 


°m,3 

c N 

c N,f 

c n> c 0 


third-moment coefficient of equivalent airfoil section 

/ x 

(see ref. 1), — =A- 

qSc 4 

resultant-force coefficient on arbitrary section, — 

qS 

resultant-force coefficient of flat plate, ^ 

coefficients of sine-series expansion of vorticity distri- 
bution on equivalent airfoil section at arbitrary depth 


c 

D 

°1 


chord 

total drag force 

drag force excluding crossflow drag and skin friction 
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D f 

d 

E 

F 

g 

L 

Ll 

l 

M 

m 3 

m = C^/a 

P 

Pc 

Po 

Pv 

q 


s 

s 

u 


drag force due to crossflow (normal to section) 
drag force due to skin friction 
leading-edge depth of submersion 

Jones ' edge correction, ratio of semi perimeter to span 
(see ref. Ik) 

resultant force 

acceleration due to gravity 

total lift force 

lift force due to crossflow 

lift force exclusive of crossflow, L - 1^ 

distance from hydrofoil leading edge to stagnation point 
in chords 

moment about leading edge 

third moment about leading edge, 2 / p(x)x^dx 

J O 


pressure, lb/sq ft 

pressure within cavity, lb/sq ft 

pressure at mean depth of hydrofoil, lb/sq ft 

fluid vapor pressure, lb/sq ft 

free-stream dynamic pressure, 

distance from image vortex to point on equivalent airfoil 
in chords 

area, sq ft 

span, ft 

perturbation velocity in X-direction 
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V 


V 

X,Y 

X 


Yl 

Yu 


a 


c 


a i 

oo 


a 


1 


r 

7 

& 

e 

9 


P 


a 


<*i 


speed of advance, fps 
perturbation velocity in Y-direction 
coordinate axes 


distance from leading edge along X-axis 

i 

distance from leading edge to center of pressure 
distance from the X-axis to hydrofoil lower surface 
distance from the X-axis to hydrofoil upper surface 


geometric angle of attack, radians unless otherwise specified 

angle-of -attack increase due to camber, ^Cq' + C^_ - C 2 / 2 ^ 
or at infinite depth a c ^ = ^A-^ - A 2 / 2 ) , radians unless 


otherwise specified 


induced angle of attack, radians unless otherwise specified 

angle between hydrofoil chord line and reference line, posi- 
tive when chord line is below reference line, radians 
unless otherwise specified 

angle of attack measured from hydrofoil chord line, 
a' = a + Oq, radians unless otherwise specified 

circulation, strength of single vortex 

central angle subtending chord of circular-arc hydrofoil 

spray thickness at infinite distance downstream 

deviation of resultant-force vector from normal to hydrofoil 
reference line 

parameter defining distance along airfoil chord, 

£ = £(l - cos 9) 

2 

lb— sec^ 

mass density, ; 

ft 14 - 

Do — Pc 

cavitation number, 

q. 

cavitation number at inception 
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0 

♦ 


n 


correction factor for variation from elliptical plan form 
angle between spray and horizontal 

angle between x-axis and line joining image vortex with a 
point on equivalent airfoil 

vorticity 

indicates "function of" for example, C^jaj = Cjj^pja + a, 


Z = X + iY 


3 


tan' 


-1 £d 

c l 


a + e 


Subscripts : 

e effective 

0 zero depth of submersion 

t total 

00 infinite depth of submersion 

c.p. center of pressure 

Barred symbols refer to equivalent airfoil section and unbarred 
symbols refer to the supercavitating hydrofoil section. 


DESCRIPTION OF SUPERCAVITATING FLOW 


The parameter defining cavity flow is 


Pn - Pr 


a = 


where p n is 


q "O 

the pressure at the mean depth, p c the pressure within the cavity, and 
q, the dynamic pressure. The magnitude of a for the condition at which 
cavitation is incipient is defined by the particular value a. . If cr 


is reduced below o±, cavitation becomes more severe; that is, the cavita- 
tion zone extends over a larger area. When a hydrofoil operates at suf- 
ficiently low values of a, the cavity formed may completely enclose the 
upper or suction surface and extend several chords downstream as shown in 
figure 1(a) . Theoretically if the cavitation parameter is reduced to zero, 
the cavity formed will extend to infinity. The flow regime where the 
cavity length exceeds the chord is defined as supercavitating flow. 
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Supercavitating flow, may be obtained by either increasing velocity 
or cavity pressure or both. At a constant depth and water temperature, 
o and therefore, the length of the cavity is dependent only on the 
velocity since p Q - p c is then p Q - p v and is constant. If part or 
all of the boundary layer of a configuration is separated, the eddying 
fluid in the separated region can be replaced by a continuous flow of 
lighter fluid such as air. (See refs. 5 and- ^«) Regulation of the amount 
of air supplied will control the cavity pressure and thus the length of 
the cavity formed. If the quantity of air supplied is very large, the 
cavity pressure will approach the ambient pressure p Q and a very long 
cavity will result even at low-stream velocities. 

The ventilation of surface piercing hydrofoils is therefore a super- 
cavitating flow due to large quantities of air supplied from the atmos- 
phere to separated flow on the suction surface of the hydrofoil. Super- 
cavitating flow as a result of ventilation also occurs when a nonsurface- 
piercing hydrofoil of moderate aspect ratio operates near* the free surface 
(see fig. 1(b)). As pointed out in reference 5> air is entrained in the 
trailing vortices and drawn to the suction side of the hydrofoil causing 
a long trailing cavity to enclose completely the upper surface of the 
hydrofoil and extend far downstream. The ventilated type cavity described 
in reference 5 differs in shape from those formed in deeply submerged flow 
because of the proximity of the free water surface. It is similar to 
planing, with the spray forming the upper surface of the cavity. Since 
the cavity pressure is approximately the same as the ambient pressure 
(at small depths of submersion), the cavitation number for this type 
flow is nearly zero. 

The experimental portion of the present investigation utilized both 
of the methods of obtaining low values of a discussed in the foregoing 
paragraphs; that is, varying velocity and cavity pressure. However, 
the case of zero cavitation number was of particular interest since it 
represents the condition of minimum lift and also because a major portion 
of the existing cavity flow theory assumes this condition. Therefore, 
most of this report treats the ventilated type supercavitating flow. 


ZERO-CAVITATION-NUMBER THEORY 


TWO-DIMENSIONAL THEORY 


Flat Plate 

The characteristics of a two-dimensional inclined flat plate in an 
infinite fluid, operating at zero cavitation number, have been obtained 
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by Kirchoff and Rayleigh (ref. 6) . The resultant force on the plate is 
given by the well-known equation 


CN,f 


2it sin a 
4 + it sin a 


( 1 ) 


The distance from the leading edge to the center of pressure is 


(X) 

= Ms 

> _ 3 

cos a \ 

Wc.p. 

4 \‘ 

4 + 

it sin a/ 


( 2 ) 


Similar work was performed by Green (refs. 2 and 7 ) to include the 
effect of the free water surface (but neglecting gravity). The solution 
is necessarily obtained in terms of the spray thickness 6 rather than 
the more useful depth of submersion and is given as two parametric equa- 
tions in terms of the parameter b. 


Cl = Cfljf cos a 



- 1/sin a cos a 
D 


ma 


(3a) 


o D - cos a /-.n \ 

c " D (5b) 

where 


2 cos a + (b cos a - l) log- ^ - 1 -^ 

e b + 1 

This result is plotted as the variation of m with &/c for various 
angles of attack in figure 2. 

Since gravity is neglected in Green's solution, it is important to 
understand the effect of the absence of this force on the result obtained. 
If the Froude number based on chord V 2 /gc is large, then near the plate, 
where the streamline radii of curvature are small, the inertia forces are 
indeed large compared to the force of gravity. Thus, under these condi- 
tions, equations (3) for the force on the plate should be applicable, for 
a given s/c, in spite of the neglect of gravity in its derivation. How- 
ever, the variable 8/c is not usually known. It is desired to know the 
force on the plate for a given depth of submersion; therefore, the rela- 
tionship between the leading-edge depth of submersion and the spray 
thickness is needed. The effect of gravity becomes very important in 
determining this relationship. Because of the neglect of gravity in the 
analysis of Green, the plate is always located infinitely above the still 
water surface - a physical impossibility in a real fluid. Gravity cannot 
be neglected in determining the location of the free water surface. The 
reason for this is that several chords away from the plate the streamline 


D = 


Jb 2 - l) 


sin a + — 

it 
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radii of curvature become very large and thus the inertia forces become 
of the same or lower order of magnitude than the force of gravity. When 
gravity is present the distance of the plate above the water surface at 
infinity is at least limited to the dynamic pressure head. The effect 
of finite aspect ratio is also to place the plate finitely near the free 
water surface as will be discussed in a subsequent section. 

In summary, if the Froude number is large, the forces on an inclined 
flat plate near a free water surface can be obtained in terms of the spray 
thickness b/c from equations (3)* The relationship between the spray 
thickness and the actual leading-edge depth of submersion cannot be deter- 
mined from Grer.i's analysis, but the relationship is known to be influ- 
enced by the presence of gravity and the aspect ratio of the plate. The 
method by which the lines of constant depth (shown in fig. 2) were obtained 
is discussed in a subsequent section. 

The angle (J) that the spray makes with the direction of motion is 
also obtained from reference 2 as 


The variation of (f> (eq. (L)) with b/c (eq. (3b)) is shown in figure 3. 


Cambered surfaces can theoretically be analyzed in two dimensions by 
the method of Levi-Civita (ref. 7). However, like many conformal mapping 
problems the method is very difficult to apply to a particular configura- 
tion and only a few specific solutions have been obtained. Among these 
is the work of Rosenhead (ref. 8) and Wu (ref. 9 ) • Although the solution 
of Wu is applicable in principle to arbitrary sections, the solution was 
presented only for the circular arc. A particular advantage of Wu's 
solution is that it includes the effects of nonzero cavitation number. 

The most useful treatment of cambered surfaces is the linearized 
theory of Tulin and Burkart (ref. 1) which is readily applicable to any 
surface configuration (with positive lower surface pressures) as long 
as the angle of attack and camber are small. The principal results of 
this linearized theory are summarized below. 

The supercavitating hydrofoil problem in the Z-plane is_transformed 
into an airfoil problem in the Z-plane by the relationship Z = -Jz. By 



( 4 ) 


Cambered Sections 
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denoting properties of the equivalent airfoil with barred symbols and 
those of the hydrofoil with unbarred symbols, the following relationships 
are derived: 


CVJ 

IX 

II 

'ix 

( 5 ) 

ti(x) = u(x^) 

(6) 

! C„ = |(ao 

(7) 

^(cl)' = |(ao +^) 

(8) 


Cm = C m , 5 = ihAo + 7Ai - 7A 2 + 5 A 3 


-%) 


( 9 ) 


The coefficients are the thin-airfoil coefficients in the sine- 

series expansion of the airfoil-vorticity distribution 


ft(x) = 


/ 00 

= 2V|Ao cot | + £ Aj, 

\ n=l 


sin n9 


( 10 ) 


where 


X = -|(l - COS o/ (0 i? 0 ^ Jt) 

and can be found for a given configuration from the following equations 

(Ha) 


A n = - ^ f — d9 + a = a + AH 

^ *J 0 dx ^ 
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An " 


^ f ^ cos n9 d0 

*J 0 ** 


(lib) 


Ey using equations (7) and (8), the condition that Aq in equation (10) 
is zero, that D(x) does not contain harmonics greater than two, and 
that ft(x) is everywhere positive over the chord (that is, positive 
lower surface pressures on the hydrofoil), Tulin and Burkart (ref. l) 
obtained a low-drag family of hydrofoils. This particular family of 
hydrofoils is given by the equation 


y _ ^1 x 8/x\^/ 2 . /x\ 2 

c 2 c 5\c ) “ \c) 


( 12 ) 


The lift-drag ratio of this foil at its design incidence (that is, a = 0) 
is from equations (7) and (8) 


L = 2 5/_2_\ 
D 4^2C l J 


(13) 


This 6.25-fold improvement over the L/D of a flat plate is most encour- 
aging and was the impetus for this investigation. 

The restrictions imposed by the assumptions of the Tulin-Burkart 
theory prevents its use in the calculation of the characteristics of 
hydrofoils suitable for use as aircraft landing gear. Here, because of 
the high hydrofoil loads on necessarily thin hydrofoils the aspect ratio 
may be as low as 1 or 2. Also the hydrofoil must operate near the free 
water surface and in some instances at large angles of attack. Thus, 
the effects of these variables on the characteristics of supercavitating 
hydrofoils (particularly of cambered sections) is needed. Much of this 
information can be obtained by certain modifications to the existing two- 
dimensional theory discussed in preceding paragraphs. 


MODIFICATIONS OF TWO-DIMENSIONAL THEORY 


Nonlinear Equation for Lift at Infinite Depth 

For any configuration the reference line of .the section from which 
the angle of attack is measured can be chosen such that Aq' = 0. It is 
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convenient in this report to assume that a = 0 refers to the reference 
line which makes Aq* = 0. With this assumption, equation ( 7 ) may he 
written as 

C L = §( a + A 1 - y) = |(«- + a c ) ( lif ) 


where a c 



is the effective increase in angle of attack due to camber 
Thus, the solution for cambered hydrofoils is merely the flat 


plate linearized solution with a replaced by a + a_. This is 

exactly analogous to the influence of camber on airfoils in an infinite 
fluid where there is an effective increase in angle of attack due to the 
camber. By carrying this procedure further, and by applying it to the 
resultant force rather than the lift, the nonlinear solution of Rayleigh 
becomes applicable to arbitrary configurations simply by replacing a 
by a + a c ; that is. 


C 




N " 


2ji sin(a + a c ) 

4 + n sin(a + a c ) 


(15) 


The lift will then be 


„ ^ 2ji sin(a + a c ) 

C L = cos 0 

4 + n sin(a + a c ) 


( 16 ) 


Here, p = a + e, where e denotes the deviation of the resultant-force 
vector from the normal to the hydrofoil reference line. For large values 
of a, e is small compared with a and cos p = cos a. When a is 
very small, about 0 ° > 

tions ( 7 ) and (8) to he 


e is a maximum and can be shown by use of equa- 
A . 2 

if Aq' = 0 as assumed. For any 


Ma, - ^ 


practical low-drag hydrofoil, this value will almost always be less than 
3° for which the cosine is very nearly 1 or cosine (a + e) = cos a = 1. 
Therefore, cos 0 in equation (l6) may be replaced by cos a with little 
loss in accuracy and great gain in simplicity. Equation (l6) then becomes 


C L 


2jt sin(a + a c ) 

4 + K sin(a + a c ) 


cos a 


(17) 
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For a circular-arc hydrofoil of central angle 7, it is shown in 
appendix A that a c = (9/16)7 and that the reference line must he chosen 
at an angle 7/8 to the chord line so that Aq' = 0. The result obtained 
by substituting this value into equation (17) is compared in figure 4 
with the linear solution of Tulin and Bur kart (eq. (7)) and- the nonlinear 
solution of Wu (ref. 9) for two circular-arc profiles. The agreement of 
equation (17) with the more exact solution of Wu is good over the entire 
range of angle of attack from 0° to 90°. Similar agreement is expected 
for any configuration of small camber. 


Nonlinear Equation For Lift at Finite Depth 

The successful modification of the Rayleigh equation to include 
cambered configurations leads at once to a similar modification of the 
solution of Green. However, in this case the argument for replacing a 
by (a + a c ) is very weak unless the section coefficients which deter- 
mine a c are known as a function of the depth of submersion. 

An examination of the linearized expressions for the lift coefficient 
of arbitrary foils at infinite depth and at zero depth reveals that both 
the lift-curve slope and the increase in angle of attack due to camber 
do change with depth of submersion. At infinite depth the linearized 
expression for lift coefficient is given by equation (7). At zero depth 
the lift coefficient must be one half of the fully wetted value obtained 
from thin-airfoil theory as pointed out in reference 10; that is, 

CL,0 = *^0,h + j ( lQ ) 


where Aq h and Aq h are the thin-airfoil coefficients of the section 
in the hydrofoil plane and are given by the expressions 


^ ,h 


1 



(19a) 


a -1 >, = - r c ° s 0 <i9 (19L) 

l,h Jt J 0 dx 

For the Tulin-Burkart section at zero angle of attack these values may 
be determined as 
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A^ h = 0 . 227 A-L (20a) 

A-^ h = 1.15IA-L (20b) 


Thus, from equations ( 7 ) and. (l8) it is seen that, for a flat plate at 
small angles, the lift coefficient goes from at infinite depth to 

na at zero depth (as given by Green) whereas, for the Tulin-Burkart 
section at zero angle of attack these values are ^(l. 25 Ap) at infinite 
depth and n(0.802Ap) at zero depth. Although the flat-plate lift coef- 
ficient doubles in going from infinite to zero depth, the ratio is only 
1.28 for the cambered section. The important point to note is that the 
value of a c for the Tulin-Burkart section changes from 1.25Ap to 

0.802A-L. 

It is now desirable to determine a c for finite depths of submersion. 
This can be accomplished by modifying the linearized theory of reference 1 
to include the effects of the free water surface. 

The effect of the free water surface may be obtained by finding the 
transformation which_will map the free water surface, the hydrofoil, and 
the cavity into the X-axis of the Z-plane. The required transformation 
is of the form Z = Z + log e Z + K 2 . By using this transformation it 

is possible, in principle, to determine the forces on the hydrofoil by 
following the procedure used in reference 1. However, the solution by 
this method became very cumbersome and was abandoned when an approximate 
method was discovered. The approximate method continues with the trans- 
formation Z = -Jz used in reference 1. 

In figure 5 it may be seen that Z = -{z transforms the free water 
surface in the hydrofoil plane, where u = 0 (see fig. 5(a)), into a 
hyperbola in the third quadrant of the airfoil plane (fig. 5(b)). The 
boundary condition that must be satisfied on this hyperbola is that the 
perturbation velocity u be zero because in complex velocity problems 
of the type considered here, it is the lines of constant velocity which 
are being transformed and not lines of constant velocity potential or 
stream function. 

For the particular case of zero depth the hydrofoil problem is trans- 
formed by Z = -yjz into the fourth quadrant of the airfoil plane. Thus, 
it may be seen in figure 5 ( 8 ) that the free water surface adds the con- 
dition that u = 0 along the negative Y-axis. This additional boundary 
condition can be satisfied (along with the other infinite-depth boundary 
conditions) by locating an image of the airf oil- vorti city distribution 
along the negative X-axis. The direction of this vorticity must be 
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opposite to that of the airfoil in order to make u = 0 at all points 
along the negative Y-axis. A simpler and often used approximation is to 
replace the distributed image vorticity by a single vortex, equal in 
strength to the airfoil circulation, at a location equal to the distance 
from the leading edge to the airfoil center of pressure as shown in 
figure 5("b) • 


For finite depth of submersion, the condition that u = 0 must be 
satisfied at all points on the hyperbola and on the negative X-axis. It 
is not possible to satisfy these conditions with a single vortex as was 
done for the case of zero depth. However, the influence on the airfoil 
of the infinite array of vortices needed to satisfy the boundary condi- 
tions shown in figure 6 may be approximated by a single vortex of strength 
T in the location shown. The adequacy of the approximation can be deter- 
mined by calculating the effect of this image vortex on the lift of a 
flat plate as the depth of submersion is varied and then, comparing the 
result with the exact solution of Green. This depth effect may be deter- 
mined by concentrating the airfoil circulation at its center of pressure 

and the image circulation at a point -■ + cJ— forward of the leading 

- rv 4 1C 

edge and £J— below the leading edge (see fig. 6) and computing the 

total downwash on the flat plate at its 22. point (see ref. ll) . The 

method assumes that the center of pressure of the airfoil remains con- 
stant at c/4 as the depth changes. The resulting downwash angle a 

at the 22. point is calculated to be 


a = 


r 

itcV 




1 + 


jdA 


,1 + 2 \jd/c + (5/4) (d/c) 


( 21 ) 


Since = ncaV, the ratio of hydrofoil lift at finite depth Cl to the 
lift at infinite depth is 


C L = = iK = T_ 

C L,°o ^m,oo ic L>00 r °° 

therefore 

C L _ _T. = 1 

C L,co r« ^ l/ 1 + /d/c 

2 \l + + (5/4) (d/c) 


( 22 ) 


(25) 
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Equation ( 23 ) is compared in figure 7 'with the exact solution of Green 
(see fig. 2, a = 0) , and the agreement is excellent. 


The adequacy of the method used in determining the influence of free 
water surface proximity on the lift of a flat plate justifies its use on 
cambered foils. However, for cambered foils the problem is more difficult 

35 

because a control point such as — for a flat plate is not sufficient to 

4 

determine the section coefficients defining the vorticity distribution. 
This final vorticity distribution fl(x) resulting from the camber and the 
image vortex must be determined, particularly if a knowledge of the pres- 
sure distribution and thus the drag and center of pressure is desired. 


In figure 8, 
on the equivalent 


it may be assumed that the final vorticity distribution 
airfoil is given by the equation 


ft(x) 


2V 


^0 



n=l 


sin 


n0 1 


(24) 


Then the induced velocity at x 1 due to fl(x) is, from thin-airfoil 
theory (see ref. 12), 


v 



n=1 


C n cos n9 


(25) 


The resulting total circulation due to the vorticity given by equation (24) 
is, from thin-airfoil theory, 

r = jtcV^Q + (26) 

Since the image vortex has a strength equal and opposite to P, the 
velocity induced by at a point x' on the airfoil is 


^i 


cos ilr 

2itrc 


(27) 
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where 


or 


where 


or 


Replacing x f 
Fourier series 


where 


r 




cos 




(28) 



a 



Cp 

Co + Cl - ■£ 

Co 


(29) 


by ^(l - cos 9 ), equation (28) can be expanded in a 
as 


00 

= f ( 9 ) = Bq + ^ ^ cos n 6 
n=l 


( 30 ) 
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and 


B 


n 


2 

Jt 



f(9)cos n9 d9 


By substituting equations (2 6 ) and ( 50 ) into equation (27) gives 


v(c 0 + %) / 

00 \ 
\ 


n - V v 2 Mbo + 

) B n cos n9] 

(31) 

\ 

h ) 


The equivalent airfoil slope ^ 

die 

when expanded in a cosine 

series 


A n cos n9 

(32) 

n=l 




Equating the resulting streamline slope to the equivalent airfoil slope 
gives 


ft = l + h = 

V V dx 


(33) 


The substitution of equations (25), (31), and (32) into equation (33) 
gi ves 


00 



n=l 


cos n9 + 



00 



n=l 


cos n9 


(34) 


By equating coefficients of like terms 9 the C coefficients are deter- 
mined as 


Aq(4 + 62 ) + AjBq 
4 + B ± - 2B q 


(35a) 
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Ci = 


2Al(2 - B 0 ) - 2^ 
4 + Bq - 2B 0 


(35b) 


- A n - 


(2Aq 4- 

^4- + B-^ — 2Bq 


(35c) 


If Aq = a and. A n = 0 where n ^ 1, the effect of depth of sub- 
mersion on the lift coefficient of a flat plate can be computed from the 
coefficients obtained from equations (35)- The values of — computed 

^L.oo 

by this method have been found to be in excellent agreement with the solu- 
tion of Green for a = 0. 


For the particular condition of Aq = 0 (the case of hydrofoils such 
as the Tulin-Burkart section at zero angle of attack) equations (35) 
become 


„ _ B 0 A 1 
L o _ : 

4 + Bq - 2B 0 

c _ 2 A l ( 2 - B o) 

1 4 + Bq - 2B 0 


( 36 a) 


( 36 b) 


Cn = A n - 


4 


Bn A l 
Bq - 2B 0 


( 36 c) 


The coefficients Bo and Bn as obtained from equation ( 30 ) are plotted 
in figure 9 against the distance to the center of pressure a for several 
depths of submersion. For the special condition of zero depth the B 
coefficients for the Tulin-Burkart section are found to be I. 296 , 0.772, 
and 0.22 for Bq, Bi, and B 2 , respectively. In making this computation, 
the final center-of -pressure location a is used; therefore, a is first 
given an assumed value, the B and C coefficients determined, then 
from the resulting C coefficients, a is calculated from equation (29) 
and the procedure repeated if necessary. For the Tulin-Burkart section 
a is found to be 0.42 for d/c =0. By using the final B coefficients 
the C coefficients are determined from equations ( 36 ) as 0.595 A q, 
0.646A-|_, and -0.6lAi for Cq, Cq, and C 2 , respectively. 
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From equation ( 7 ) the ratio of the Tulin-Burkart section lift at 
zero depth to the lift at infinite depth is 


C L ,q = (0-393 + 0-646 + O. 5 O 5 ) Ax 

C L,oo 5 AiA 


(37) 


The value 1.24 compares favorably with the more exact value of 1.28 given 
in the paragraph following equations (20) . The results of calculating 

Cl 

— for finite depths of submersion for the Tulin-Burkart section are 

Cl, oo 

plotted in figure 10. Similar calculations have been made for a circular 
arc, and the configurations given in reference 13 . The results of these 
calculations are also shown in figure 10 . 

The true linearized lift-curve slope m for finite depths of sub- 
mersion in the equation Cl = m(a + a c ) is that shown in figure 2 for 
a = 0. Therefore, the effective angle of attack due to camber a c is 
obtained from the following relationship: 


m (o=0) a c _ C L 


2 a c,°o 

Therefore 


a c,<» 2m( a;= o) Cl^oc, 

a c . 

Values of are plotted against d/c in figure 11 for the Tulin- 

a C,oo 

Burkart, the circular arc, and the sections given in reference 13 . 

Equation ( 38 ) is obviously limited by the linearizing assumptions 
made in its derivation. An important limitation is due to the assumption 
that the free surface is always horizontal and thus S/c = d/c. At 
small depth-chord ratios and particularly for large magnitudes of camber 
the free water surface is not horizontal and b/c / d/c. Thus, for small 
values of d/c and large magnitudes of camber the values of a c /a c 
given in figure 11 are probably too low. 

With a knowledge of the angle of attack due to camber a c at finite 
depths of submersion. Green's solution is now modified to include camber 


C L,co 


c L 
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by treating the effective angle of attack as (a + a c ), where a c is 
obtained from figure 11. This is exactly the method used in modifying 
the Rayleigh equation to obtain the nonlinear approximation for the lift 
coefficient at infinite depth. With this assumption, the resultant-force 
coefficient for a cambered hydrofoil at any positive depth of submersion 
is obtained in terms of the spray thickness b/c from equations (3) as 


C N (a} 2 C N>f ja + a c J 


(39) 


Equation (39) states that the resultant force on a cambered section is 
approximated by replacing a in Green's solution for a flat plate by 
the effective angle of attack (a + a c ) . It will be shown that the 
resultant force will deviate only slightly from the normal (as previously 
pointed out for the condition of infinite depth) and therefore 

C L {aJ - % f| a + a, c |cos a (40) 

Actually, by itself, equation (40) is of little practical value 
because Cj^f is given in equations (3) in terms of the usually unknown 
parameter b/c instead of the more useful d/c. Thus, it is necessary 
to determine the relationship between S/c and d/c or at least an 
empirical substitute. 


Relation Between Depth of Submersion and Spray Thickness 

Certain relationships between the leading-edge depth of submersion 
and the spray thickness of a flat plate can be stated and are given 
as follows: 

(1) The trivial but useful case for a = 0 where d/c = b/c 

(2) The case of a = 90°, &/c = <», where the stagnation streamline 
is parallel to the direction of motion and in order to satisfy the con- 
tinuity equation b/c = d/c + 0.5 (see fig. 12(a)) 

(3) The case of a = 90°, b/c = 0 where again the stagnation line 
(this time the free water surface) is parallel to the direction of motion 
and b/c = d/c + 1 (see fig. 12(b)) 

Cases (2) sind (3) suggest that the stagnation line for a = 90° is 
also parallel to the direction of motion for finite b/c (as shown by 
the solid line in fig. 12(c)). However, this is not consistent with the 
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results of Green's analysis where it can he shown from momentum considera- 
tions that the streamline is curved as shown by the broken line In fig- 
ure 12(c). However, when the effects of gravity and even slight devia- 
tions from infinite aspect ratio are considered, the true location of the 
stagnation line must lie between these two lines. The tendency of gravity 
to cause this change has been discussed already. Possibly more important 
than gravity is the influence of finite aspect ratio. 

If the supercavitating hydrofoil is replaced by a system of horseshoe 
vortices, it can be shown that (for finite angles of attack and infinite 
aspect ratio) the stagnation line infinitely forward of the hydrofoil is 
also infinitely below the stagnation point - as obtained by Green. How- 
ever, for an aspect ratio of as much as 100, the stagnation line approaches 
a finite asymptote about 0.2 chord below the stagnation point when the 
lift coefficient is about 0.5* At a lift coefficient of 0.5, the asymp- 
tote becomes about 0.1 chord for moderate aspect ratios. This was observed 
in the experiments described in a subsequent section where it was noted 
for the angle-of-attack range investigated that the free water surface 
was practically undisturbed only slightly forward of the hydrofoil leading 
edge. Therefore for an angle of attack of 90° where the total circulation 
is zero, it is not unreasonable to expect the solid stagnation line in 
figure 12(c) to be more nearly correct than the dashed one when both finite 
aspect ratio and gravity are present. It should be noted that although 
the stagnation line may deviate from the theory as given by Green, that 
the pressure distribution and force on the plate are assumed to remain 
unaltered from the expressions given in equations (3) . Therefore, if at 
90° the stagnation line is assumed parallel to the direction of motion, 
it follows that for a = 9O 0 


6/c = d/c + Z/c 


(41) 


where Z/c is the dimensionless distance from the leading edge to the' 
stagnation line (see fig. 12(c)). This distance can be obtained from 
Green's work and is, for a = 90° 


1 + ^ - 
2 

" io& 4 . 1 ! 

1 - yl 

7h 

1 

OJ 

Up 

it -”- 1 k) 

Jtl 

(b - \i 

^ - 1) 

1 + log e(. 

rf) 

1 


(42) 


where b has been previously defined in equations (3). 
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Case (l) and equations (4l) and (42) give the end points (at a = 0° 
and 90 °) of the lines of constant d/c which are superimposed on fig- 
ure 2. The trends of these curves at intermediate angles of attack were 
obtained from experiment and will he discussed in a subsequent section. 

The preceding discussion in this report has been concerned with the 
effects of gravity and finite aspect ratio on the relationship between the 
spray thickness and the leading-edge depth of submersion. It is also of 
interest to examine at least qualitatively the effects of gravity and 
finite aspect ratio on the spray angle given by equation (4) and plotted 
in figure 3 . In Green's analysis, it is shown that at a = 90° the total 
force on the entire plate is F = p 6 V^(l - cos (fi ) ; that is, the force is 
exactly equal to the change in momentum due only to the flow in the spray. 
This means that the stagnation line must be curved as pointed out in the 
preceding discussion. The straight stagnation line (which leads to equa- 
tion ( 4 l)) implies that the change in momentum in the spray flow results 
only in the force on that portion of the plate above the stagnation 
point. Since the change in momentum of the spray flow now represents 
only a portion of the total force on the plate, and since 5/c is con- 
sidered to be the same for both stagnation-line curvatures , the angle (j> 
must decrease from the value given by Green in figure 3 . That is, the 
angle will change so that only the force on the plate above the stag- 
nation line is given by the expression F = p 6 V^(l - cos (ft ) . 


THREE DIMENSIONAL THEORY AT FINITE DEPTH 
Lift 

The flow about a supercavitating hydrofoil may be constructed by a 
suitable combination of sources and vortices. The vortices contribute 
unsymmetrical velocity components and lift, the sources contribute 
symmetrical components which provide thickness for the cavity but no 
lift. For a finite span the vortices can not end at the tips of the 
foil and a system of horseshoe vortices must be combined with the sources 
to describe the flow. If it is assumed that the influence of finite span 
on the two dimensional lift coefficient is due to the effects of the 
trailing vorticity then the resulting effect of aspect ratio is exactly 
the same as for a fully wetted airfoil. Jones (ref. 14) gives the lift 
of a fully wetted elliptical flat plate as 

C L = | 2 * (a - 04 .) (43) 

where E is the ratio of semiperimeter to the span and 04 is the 

induced angle of attack caused by the trailing vorticity. Thus the effect 
of aspect ratio is to decrease the two dimensional lift curve slope by a 
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factor l/E and to decrease the effective angle of attack by an incre- 
ment . Therefore for the finite aspect ratio supercavitating hydrofoil 

at infinite depth, equation (l4) is modified to give 

C L,1 = | |( a + a c - a i) (^) 

or more generally for finite depth, equation (40) becomes 

C L, 1 = | c Nf( a + a c - a i) cos a (k5) 

where for rectangular plan form of aspect ratio A, E = A - + ^ and 

A 

04. = ^(1 + t ) (46) 

-*■ itA 

where t is a correction for plan form (see ref. 12). 

Another effect due to finite aspect ratio is the concept of additional 
lift due to crossflow (see refs. 15 and l6). This crossflow lift is 
assumed due to the drag on the hydrofoil contributed by the component of 
free -stream velocity normal to the hydrofoil. In the present case of zero 
cavitation number, the crossflow drag coefficient is the Rayleigh value, 
0.88. Since this lift is caused only by the spanwise flow (flow around 
the ends of the plate) it is also modified to account for the aspect ratio 
by the Jones' edge correction, l/E. Since only the spanwise flow is con- 
sidered, E is now the ratio of semiperimeter to chord. Because the flow 
being considered is normal to the plate, the induced angle for this flow is 
zero. Thus for a flat plate, the crossflow lift C T „ is 

i-J y C 

= a + ^ 0.88 sin^ a cos a (4 7 ) 

No experimental or theoretical information on the crossflow lift of 
cambered surfaces is available in the literature. In order to approximate 
this component the following assumptions are made: 

(1) The crossflow force acts normal to the hydrofoil chord line. 

(2) The effective direction of the free stream on the plate is altered 
by the increase in angle of attack due to camber a c . 

Thus, the crossflow lift on cambered sections is assumed to be 
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C T = — — — 0.88 sin 2 (a' + a c ")cos a' (18) 

L > c A + 1 v 1 

where a' = a + a, 0 , and a 0 is the inclination of the chord line to the 
reference line of the section (positive if the chord line is below the 
reference line) and a c is obtained from figure 11 for the depth of 
interest. 

The total lift on a finite aspect ratio hydrofoil operating near the 
free water surface - f s therefore obtained by adding equation (18) to equa- 
tion (I5) "to give 



~ A 
A + 1 


C n ,f l a + a c 



cos a + 


1 

A + 1 


0.88 sin 2 (a' 


+ a c) 


cos a 


I 


(49) 


In view of the very approximate nature of equation (18) it is desir- 
able to examine the effect of this crossflow term on the total-lift coef- 
ficient. For the Tulin-Burkart , aspect-ratio-1 section (Aq = 0.2) 
used in this investigation, at d/c = O.O7I the ratio of the calcu- 
lated crossflow lift, Cl jC to the calculated total lift was 0.157 

at a = 1° and O.283 at a = 20°. For the five-term section of refer- 
ence 13, where Aq is 0.075, the aspect ratio is 3, and d/c is 0.071, 
the ratio has been calculated as 0.011 at a = 1° and 0.072 at a = 20°. 
Thus any inaccuracies in the crossflow lift as computed by equation (18) 
will appreciably affect the total- lift coefficient at large angles, small 
aspect ratios, and large cambers. On the other hand at higher aspect 
ratios and small cambers, errors in the crossflow component do not greatly 
influence the total lift. 

Equation (I9) may be written in terms of the slope m (given in 
fig. 2) as 


c L = 


A + 1 


m(a 


+ a r 


aq) 


cos a 


cos (a + a c - aq) 


— - — 0.88 sin 2 (a' + a c )cos a' (50) 

A + 1 c 

where a c is obtained from figure 11 for the depth-chord ratio of 
interest and aq is obtained from equation (16). In equation (16), 

C L q is the first term in equation (50). Equation (50) is solved by 
iteration and the convergence is quite rapid. A sample calculation of 
the lift coefficient of a cambered hydrofoil is given in appendix B. 
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Drag 

The drag coefficient of a supercavitating hydrofoil of finite aspect 
ratio operating at zero cavitation number and finite depth of submersion 
is 


Cd = Cl^I tan (a + O + Cp^ c tan a* + Cf (51) 

where is the first term in equation ( 50 ) and e is the deviation 

of the resultant-force vector from the normal. For a flat plate e = 0, 
a f = a, and thus Cp = Cp tan a + Cf . For cambered surfaces similar to 
the circular-arc or Tulin-Burkart section, e becomes very small at large 
angles of attack and may be neglected; however, at small angles of attack, 
the effect of e on the drag coefficient cannot be neglected. An approxi- 
mation to the value of e can be made by determining its value from the 
two-dimensional linearized solution and then modifying the result for the 
case of finite angles of attack and aspect ratio. 

The linearized drag coefficient as given in reference 1 is 

C D (52) 

J 0 c v 

replacing x by 4^-(l - cos 9) and dx by sin 9 d9 equation (52) 
becomes 


C D = “ 



(l - cos 9) sin 9 uv 
V 2 


d9 


(53) 


Now u may be written in terms of the vorticity on the equivalent air- 
foil operating at finite depth as 


u 




+ cos 9\ 
sin 9 / 


C n sin n9 


n=l 


(54) 
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and 


v = V 


'*0 + I 


cos n0' 


n=l 


(55) 


Therefore 


C D 


■ -/ v 

J 0 


cos 9) sin 9 


1 

c °{Hnnr) + I Cn sln ne 


cos n0 d0 


(56) 


For the condition where Aq = a. Op becomes after integrating 


Cn = 5a C n + Ct - -A + 5 


Cl £2 
k 2 


- ^)A X + 


C 0 + C 1 + A k C 


L u 2 


(57) 


At infinite depth (Cq = Aq, C n = A n ) equation (57) reduces to the 
value given by Tulin-Burkart in reference 1. 

- >2 


v-ifo 4) 

For Aq = cl, that is Aq 1 = 0 

§*(c 0 + Cl - + § _(tt - t) a i + 


C D 

a + e = — ^ = 

Cl 


(58) 


Co + Cl + -tM AkC? 

■ 2 - -h - ttj 




(59) 


Therefore 


(Cl - 2 C 2 )ai + 

( 2 C 0 + 2Ci + CiJa 2 - A^C 2 

,(< 

^0 + Ci - 

1 


( 60 ) 


The value of e given by equation (60) is adequate only for the case of 
small angle of attack and camber and depth-chord ratios larger than 
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about 1. At small spray thickness to chord ratios it may he seen in 
figure 13 that the spray angle becomes quite large even for small angles 
of attack. When such a flow is transformed by Z = -^Z" the cavity 
streamline and the free water surface are rotated as shown in figure 13. 
In figure 13(b) the boundary conditions are now different from the simple 
u = 0 used in the small-angle theory. If a system of vortices could be 
located to satisfy the boundary conditions along these new lines, a solu- 
tion for the resulting vorticity on the foil could be obtained. Such a 
method would involve taking a different spray angle (from fig. 3) for 
each depth and angle of attack in order to locate the image vortex or 
array of vortices. Also, for large angles of attack and spray angles 
the linearizing assumption of v « V is not adequate. In calculating 
the lift, these difficulties were avoided by using Green's solution which 
takes the effect of the spray angle into account. It can be seen in 

c 

figure 2 that as the angle of attack increases, the ratio of — dimin- 

Cl,o° 

ishes, therefore, the image vortex must have less influence on the 
resulting hydrofoil vorticity for large angles of attack. An approxima- 
tion to the correct hydrofoil-vorticity distribution for finite angles 
of attack operating near the free water surface can be obtained by using 
the model shown in figure 8, and increasing the value of d/c used so 

c L 


that the resulting 


' , L,oo 


corresponds to that given by Green. Thus, for 


large angles of attack, camber, and finite aspect ratio an effective 
depth of submersion, (d/c) e should be used in equation (60) . The value 
of (d/c) e is the value of d/c on the a = 0 line in figure 2 corre- 
sponding to the value of m = mg, where 


L £ 

Cl, 00 2 


m{a + q c - Oj} „ 

m {q + <*C - Oil* 2 


(61) 


The value of the C coefficients are then determined for (d/c) e and 
Aq = q - q^ 

A sample calculation of the drag coefficient is given in appendix B. 

It has been found after several calculations that the value of e and 
thus the drag coefficient is not greatly affected by the depth of sub- 
mersion. In fact, a rough approximation to e may be obtained by assuming 
C 0 = q - q ± , and C n = in equation (60) . 
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Center of Pressure 

The linearized expression for the center of pressure of a finite- 
aspect-ratio, supercavitating hydrofoil operating at zero cavitation num- 
ber and finite depth of submersion is 

_ c m,l 
x c.p.,l C L)1 


where the C coefficients are determined at the effective depth of sub- 
mersion given by equation (6l) and for A q = a - oq. Superimposed on 
this flow is the crossflow component of lift which is assumed to be 
distributed ■uniformly over the chord and acting in a direction normal to 
the chord line. Thus, the distance from the leading edge to the center 
of pressure of the crossflow -lift component *c.p.,c is given by 


JL_ 

16 


5 Cq + 7 Cq - 7C2 + 3C3 - -g- 


C 0 + Ci - — 


^2 

2 


(6 2 ) 


x c . p . , c 


0.5c 


(63) 


Admittedly, this assumption is crude and accurate only for a flat plate. 
For cambered surfaces the crossflow will not be uniformly distributed 
and for low drag cambered sections is probably concentrated on the rear- 
ward portion of the hydrofoil. 

By combining equations (62) and (63) the center of pressure of the 
combined flows is, therefore, 


x c.p. 


^L, l x c.p. ,1 + 

l C L 


0.5C L ,c\ 

— r 


(64) 


A sample calculation of the center of pressure is given in appendix B. 
As in the case of e, a few calculations reveal that a fair approximation 
for x c<p> ^i is obtained by using Cq = a - aq and C n = A n in equa- 
tion (62) . 


EXPERIMENTAL INVESTIGATION 


MODEIS 


Two models of 7.071-inch chord were used in the experimental inves- 
tigation. As shown in figure 14 (a), the first had a lower surface 
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conforming to the Tulin-Burkart low-drag configuration given by equa- 
tion (12) . The two-dimensional design-lift coefficient was selected as 
0.392 corresponding to a value of 0.2 for the coefficient A]_. Since the 
hydrofoil is designed to operate in a cavity, the shape of the upper sur- 
face is arbitrary as long as it does not interfere with the formation of 
the cavity from the leading edge. Since the greatest advantage is to be 
obtained at small angles of attack and thus thin cavities, the thickness 
of the hydrofoil must be small. For the present cambered model the upper 
surface profile from leading edge to midchord was arbitrarily chosen to 
conform with the free streamline leaving the leading edge of a two- 
dimensional flat plate at 5° incidence. (See ref. 1.) The thickness of 
the portion rearward of the midchord was made in the image of the forward 
portion and resulted in a symmetrical thickness distribution with a maxi- 
mum thickness -chord ratio of 3.3 percent. Since the center of pressure 
of the Tulin-Burkart hydrofoil is located near the midchord, this symmet- 
rical section minimizes the torsional moment on the foil and results in 
less twist than would be experienced by a nonsymmetrical section. 

The second model was of triangular cross section with a flat bottom 
as shown in figure lL(b) . The maximum thickness was 5 percent of the 
chord. 

Such thin sections lead to structural limitations in the aspect ratio 
when supported by a single strut at midspan. Since an aspect ratio of 1 
is about the most desirable from the structural standpoint and it also 
represents the accepted dividing line between hydrofoils and hydroskis, 
both models were made with a square plan form. 

The strut, which can also be seen in figure lL, had an NACA 66q-012 
airfoil section. The strut was mounted perpendicular to the flat plate 
and perpendicular to the X-axis of the cambered surface. The intersection 
of the strut and upper surface was without fillets. Both the hydrofoil 
and the strut were made of stainless steel and were polished to a smooth 
finish. 


APPARATUS AND PROCEDURE 


Tests were made by using the carriage and existing strain-gage bal- 
ances in the Langley tank no. 2 which independently measure the lift, 
drag, and pitching moment. Figure 15 shows a view of the test setup with 
the cambered hydrofoil and the balance attached to the carriage. The 
moment was measured about an arbitrary point above the model and the data 
thus obtained were used to calculate the moments about the leading edge. 
The positive directions of forces, angles, and moments used in presenting 
the force data are shown in figure 16 . 
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The force and moment measurements were made at constant speeds for 
fixed angles of attack and depths of submersion. The depth of submersion 
is defined as the distance from the undisturbed water surface to the 
leading edge of the model. 

Ventilated or zero-cavitation-number tests were made on both models 
at a depth of submersion of 0.5 inch throughout a range of angle of attack 
from 6° to 20° for the flat plate, and 8° to 20° for the cambered hydro- 
foil. Three methods of obtaining ventilated flow at this 0.5-inch depth 
were used: (l) Normal ventilation through the trailing vortices as 

described in reference A, (2) injection of air through the port on the 
strut leading edge (see fig. 14) (this air was supplied at a rate of 
0.012 pounds per second and cut off after ventilation was established), 
and (5) a l/32-inch-diameter wire was soldered on the leading edge of 
the cambered model to cause local separation and thus ventilation. Both 
models were also investigated at a depth of submersion of 0 inch at 
a = A 0 for the flat plate, and a = 6°, 8°, and 10° for the cambered 
surface. At angles of attack of 16° and 20° for the flat surface, forces 
were also measured over a range of depth of submersion for which venti- 
lation could be obtained (d = 0 to 2 inches) . 

The thickness and direction of the jet or spray leaving the leading 
edge of the flat plate were also measured at 16° and 20° for a range of 
d from 0 to 2 inches. A schematic drawing of the instrument used for 
measuring the spray thickness and direction is shown in figure 17. The 
stagnation tube was lowered through the spray during a test and the pres- 
sure and location of the tube center line recorded on an oscillograph. 
Almost instantaneous response of the stagnation-tube — pressure-cell 
combination was obtained by completely filling the tube and connecting 
line with water. The point of entering and leaving the spray was obtained 
by comparing the location of the tube with the rise and fall of pressure 
as the tube passed through the spray. The vertical location of the tube 
was obtained from the output of the slide-wire circuit also shown in fig- 
ure 8. During a test the tube was passed through the spray several times 
and an average of the results was taken. 

Tests of the nonzero-cavitation-number case were made at a depth of 
submersion of 6 inches where vortex ventilation did not occur. Measure- 
ments of lift, drag, and moment were obtained for a range of velocities 
from 20 to 80 feet per second at angles from 16° to 20° for which long 
trailing air or vapor cavities could be obtained. Data for air-filled 
cavities were obtained by introducing air from an external metered supply 
to the upper surface through the ports on the strut leading edge shown 
in figure 14. During these tests the pressure in the cavity formed, was 
measured by a pressure cell connected to a l/l6-inch-diameter orifice 
located near the bottom of the strut so as to be within the cavity. (See 
fig. 14.) This measured pressure was used in computing the cavitation 
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number for the cavity formed. The airflow rate was measured by an 
orifice-type flow meter. Airflow rates up to 0.012 pound per second 
were obtained with the test arrangement. 


ACCURACY 


The change in angle of attack, due to structural deflection caused 
by the forces on the model was obtained during the calibration of the 
balances and test data were adjusted accordingly. The maximum correction 
necessary was only 0.1°. The estimated accuracy of the measurements is 
as follows: 


Angle of attack, deg ±0.1 

Depth of submersion, in ±0.1 

Speed, fps ±0.2 

Lift, lb ±15 

Drag, lb ±7 

Moment, ft-lb ±6 

Cavity pressure, lb/sq ft ±10 

Spray thickness, in ±0.05 

Spray angle , deg ± 1.5 


The forces and moments were converted to the usual aerodynamic coef- 
ficient form by using a measured value of the density of 1.93 slugs per 
cubic foot. The kinematic viscosity measured during the tests was 

1.70 X 10 - ^ pounds-second per square foot. Thus, for the range of veloc- 
ities investigated, the Reynolds number based on chord ranged from 

0.7 x 10 6 to 2.8 x 10 6 . 


RESULTS AND DISCUSSION 


ZERO CAVITATION NUMBER 


General 

Ventilation inception .- The process of ventilation by air entrainment 
through the tip vortices of a hydrofoil near the free water surface is 
described in reference 4 and is shown for the present models operating at 
a depth of 0.5 inch, in figures 18 and 19 . The speed at which complete 
ventilation occurred, when tested at this depth is shown in figure 20 . 
Complete ventilation by entrainment through the vortices was not possible 
at angles less than 10 ° for the flat surface and 12 ° for the cambered 
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surface. Figure 18(a) shows the incomplete ventilation of the flat plate 
at 8 ° and figure 19 (a) the similar action of the cambered surface at 10 °. 
Figure 18(h) shows the complete ventilation of the flat surface at 10° and 
figures 19 (b) and (c), the complete ventilation of the cambered surface 
at 12° and l 6 °, respectively. Leading-edge cavitation of both models may 
be noted in the photographs. 

One of the requirements for establishing ventilated flow from the 
leading edge is that the hydrofoil upper surface lie beneath the upper 
cavity streamline f~r the condition of zero cavitation number. Thus, 
for a given hydrofoil thickness there is a minimum angle below which 
ventilation from the leading edge cannot exist. Figure 21 shows the 
theoretical infinite depth location of the upper cavity boundary for two- 
dimensional hydrofoils having the same cross section as the models tested. 
These streamlines, computed from equations given in references 1 and 1J, 
reveal that ventilation of the two-dimensional flat plate should be 
possible at angles greater than 3 * 2 °. 

As pointed out in references 1 and 13, the slope of the cavity 
streamline is determined by the magnitude and distribution of the circu- 
lation on the equivalent airfoil section. The reduction of the two- 
dimensional (infinite depth) circulation by an amount equivalent to 
Jones' edge correction, l/E = A/(A + l) , will cause the cavity streamline 
to be lowered over the entire span by an amount proportional to l/E. 

When the induced angle <x^ is considered, the streamline is brought even 
closer to the hydrofoil upper surface, particularly near the tips. 

By considering these two effects and by assuming that the influence 
of ai_ is uniform over J -he span, the effect of finite aspect ratio on 


Thus, in order to establish a cavity at angle of attack a and aspect 
ratio A which is equivalent to the cavity formed at angle of attack a 
and A = 00 the following equation must be satisfied: 


the cavity slope 



may be approximated by the following equation: 



( 65 ) 


a 



( 66 ) 
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By using equation (66) and the value of (a) = 3*2°, the geometric 
angle of attack of the aspect-ratio-1 plate would have to he about 8° 
in order to ventilate from the leading edge at infinite depth of sub- 
mersion. The effect of the proximity of the free water surface is to 
rotate the cavity streamline away from the hydrofoil as seen in fig- 
ure 3. Thus, it should he possible to establish a vented cavity on the 
flat section at angles at least as small as 8° when operating near the 
free water surface. Therefore, the absence of ventilation at angles 
less than 10° cannot be attributed to interference from the model upper 
surface. In fact the formation of leading-edge cavitation at 8° indi- 
cates that if the velocity could have been increased well above 80 feet 
per second (preferably to zero cavitation number) a vapor cavity com- 
pletely inclosing the upper surface would result. The reason that venti- 
lation did not occur at angles less than 10° could be that the region of 
boundary-layer separation in the midportion of the chord did not extend 
to the aerated vortices. This conjecture was proved to be the case by 
injecting air from the external supply down the strut to the upper sur- 
face of the model. Complete ventilation of the upper surface was then 
possible at angles as low as 6°. The air was supplied (at the maximum 
rate* of 0.012 pound per second) during the acceleration of the model 
and shut off after ventilation was established. The process is shown 
in figure 18(c) for the flat plate at 8° angle of attack. 

Figure 21 shows that ventilation of the two-dimensional cambered 
section should be possible at angles above about k°. In reference 13 
at large angles of attack, the cavity shape is almost independent of the 
camber, depending only on the angle of attack. Therefore, the finite- 
aspect-ratio corrections given by equations (65) and (66) are assumed 
applicable to cambered sections. When these corrections for finite 
aspect ratio are taken into account the minimum angle for ventilation of 
the cambered section becomes about 12°. Again the proximity of the free 
water surface will reduce this infinite-depth estimate. Thus, failure 
to establish complete ventilation of this hydrofoil at angles less than 
12° (at d = 0.5 inch) may have been due to interference of the upper 
boundary of the model. In fact it was founr? that injecting air to the 
upper surface at angles less than 12° failed to establish ventilation 
or to influence the ventilation speed at 12°. Since the lift-drag ratio 
of cambered supercavitating hydrofoils inc^ oases with decrease in angle 
of attack, these observations suggest the use of higher geometric aspect 
ratios or the effective increase of aspect ratio by the use of end plates. 
End plates appear to be particularly applicable if low-aspect-ratio 
cambered hydrofoils are to be designed to ventilate from the leading edge 
at small angles of attack and large depths of submersion. 

Since complete ventilation seemed to require extensive boundary- 
layer separation on the upper surface, a l/32-inch diameter wire was 
soldered on the leading edge of the cambered model to cause a local 
boundary- layer separation. With this wire in place it was possible at 
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the 0.5-inch depth of submersion to obtain complete vortex ventilation 
of the upper surface at 10° incidence but not at 8° incidence even at the 
maximum available velocity of 80 feet per second. However, by introducing 
air through the ports on the strut it was possible to obtain complete 
ventilation at angles as low as 8°. Evidently air brought up by the tip 
vortices at 8° incidence could not reach the small separated region behind 
the wire; however, the forced air was capable of reaching this zone and 
establishing a cavity. The cavity formed extended downstream far enough 
to intercept the ventilated tip vortices and complete ventilation then 
occurred. Of course, the primary influence of the wire was to increase 
the cavity ordinates so that the hydrofoil upper surface did not inter- 
fere with establishing ventilation from the wire. Although the basic 
force data obtained with the wire were recorded and are presented, they 
axe not used in the later correlation with theory because of the unknown 
influence of the wire on the forces. 

Basic force and moment results .- The basic data from the tests of 
the ventilated hydrofoils at a depth of submersion of 0.5 inch are pre- 
sented in figure 22 for the flat plate and figure 23 for the cambered 
surface as curves of lift, drag, and pitching moment about the leading 
edge against speed, for various angles of attack. Ventilated-flow data 
obtained at zero depth of submersion at L° incidence for the flat plate 
and 6°, 8°, and 10° for the cambered surface are also included. 

The basic ventilated flow data obtained for depths of 0, 1.0, 1.5* 
and 2.0 inches for incidences of 16 ° and 20° are presented in figure 2b. 


Comparison of Experimental Results With Theory 

Spray thickness . - In the section on theory the need for determining 
the relationship between the leading-edge depth of submersion and the 
spray thickness for a flat plate was pointed out. These variables were 
measured for the flat plate at 16 ° and 20° over the range of depth of 
submersion from 0 to 2 inches. The data are shown in figure 25 . The 
spray-thickness measurements presented were obtained approximately 20 per- 
cent of the chord rearward of the leading edge. On one test, measurement 
of the spray thickness was also obtained at about the midchord and the 
results were in agreement with those obtained at the 20-percent location. 
Therefore, for the range of depths investigated, the spray thickness 
measured can be considered as the theoretical value infinitely rearward 
of the foil. It may be noted in figure 25 that the spray thickness is 
greater than the leading-edge depth of submersion, the magnitude of the 
ratio increases with either increase in angle of attack or decrease in 
depth of submersion. Also shown in figure 25 are the theoretical rela- 
tionships between d/c and d/c for the two-dimensional 0° and 90 ° cases 
previously discussed. It may be noted that the trends of the experimental 
and theoretical curves are parallel. 
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In foregoing sections it has been pointed out that the relationship 
between leading-edge depth of submersion and spray thickness is influ- 
enced by the force of gravity and thus cannot be obtained from Green's 
analysis. Also the influence of aspect ratio has been shown to be impor- 
tant. The relationship for a = 0° and 9O 0 has been determined but the 
manner in which intermediate angles of attack affect the result was not 
understood. The data shown in figure 25 combined with the end-point 
results previously established for a = 0° and 9O 0 permit lines of con- 
stant d/c, at least for the aspect-ratio-1 condition, to be drawn on 
figure 2. The experimental data shown in figure 25 were plotted on fig- 
ure 2 at their equivalent angle of attack a - cuj_. This corresponds to 
equivalent angles of attack of approximately 12. 5 0 for the 1 6° case and 
15. 9° for the 20° case for the range of d/c presented. For values of 
d/c greater than O.285, only the theoretical end points were available 
and paralleling lines were faired in. Although these lines are accurate 
only for the aspect-ratio-1 condition, they are considered to be good 
approximations even for aspect ratios as high as 6. Such an approxima- 
tion is reasonable because the asymptotic value of the stagnation line 
infinitely forward of a flat plate is only of the order of about 0.1 chord 
below the stagnation point for aspect ratios less than 6 and lift coef- 
ficients less than 0.5- It is obvious from figure 2 that for depths 
greater than 1 chord, the end points may be connected by any reasonable 
line (for example, a straight line) with very little loss in accuracy. 

Spray angle .- Figure 26 shows the effect of depth of submersion on 
the spray angle 0 for the flat plate at l6° and 20°. The theoretical 
spray-angle variation also shown in figure 26 was obtained from figure 3 
by considering the two-dimensional angle of attack as a - ai and the 
two-dimensional d/c as the actual measured value. The tangent of the 
spray angles obtained from figure 3 for the angle a - were then 
reduced by the factor l/E. The justification for this l/E modification 
is the same as that presented in the section on ventilation inception. 

That is, that the slope of the upper cavity streamline is proportional 
to the equivalent airfoil circulation and thus if the two-dimensional 
circulation is reduced by l/E, the two-dimensional spray-angle slope 
will also be reduced by this amount. When these three-dimensional cor- 
rections are made, the calculated angles are still high at small values 
of d/c and appear to coincide with the experimental data for values 
of d/c greater than about 0.25- Since the theorem that the slope of 
the cavity streamline is directly proportional to the circulation was 
obtained in reference 1 for very thin cavities at infinite depth; it 
cannot be expected to be applicable when the spray angles are very large. 
This deviation from assumed small angles may be the only explanation 
required for the discrepancy between theory and experiment at small values 
of d/c. However, it might also be attributed to the influence of gravity 
and finite aspect ratio on Green's solution. In the section entitled 
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"Relation Between Depth of Submersion and Spray Thickness," it was sug- 
gested that the spray angles given in figure 3 would be reduced because 
of these effects and at least this assumption has not been disproved by 
the experimental data. 

Lift coefficient .- All ventilated force and moment data in coefficient 
form were found to be independent of speed in the range tested. (This 
independence was true in the present investigation because of the shallow 
depth of submersion and therefore a = 0. At very large depths, a will 
be greater than zero because p Q is greater than p c even if the cavity 
is fully vented to the atmosphere. Therefore at large depths, changes in 
velocity will affect the lift coefficient because these changes affect cr.) 
The data shown in figures 22(a) and 23(a) are plotted in figure 27 as lift 
coefficient against angle of attack for each of the models tested. A 
comparison of the lift coefficients of the two models shows an effective 
increase in angle of attack of the cambered model as predicted by the 
Tulin- Burkart theory. Also shown in figure 27 are the theoretical lift- 
coefficient curves obtained from equation (50) . The theory is about 
3 percent lower than the measured values. Thus, the use of equation (50) 
in engineering calculations of the lift coefficient appears to be warranted. 

The variation of lift coefficient with depth of submersion of the 
flat-plate model at angles of incidence of 16° and 20° is shown in fig- 
ure 28. Note the slight increase in lift coefficient as the hydrofoil 
nears the surface. Also shown in figure 28 is the theoretical variation 
of the lift with depth of submersion obtained from equation ( 50) • The 
theory is in excellent agreement with the data and accurately predicts 
the increase in lift with decreasing depth of submersion. 

Further verification of equation (50) is shown in figure 29 where a 
comparison is made with the experimental data of Fuller (ref. 17) . These 
data were obtained on sections with lower surface profiles similar to 
those of the present investigation but at an aspect ratio of 2. The 
calculated values are in good agreement with the measured values for 
both the flat and cambered models. 

Drag c oef f i cient . - In figure 30 the data of figures 22(b) and 23(b) 
in coefficient form are compared with theoretical values obtained from 
equation (51). The friction drag coefficient of one side of either of 
the models was calculated to be about 0.003- With the strut drag included, 
the total-skin-friction drag coefficient Cf was taken as 0.004. By 
using this value of Cf in equation (51), the agreement between theory 
and experiment is good for both models. 

In figure 31 the experimental lift-drag ratios obtained from the 
data of figures 22 and 23 are compared with theory. Again, both experi- 
ment and theory include the skin-friction drag coefficient- (Cf = 0.004). 
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The agreement between theory and experiment is good for both models. 

The superiority of the cambered hydrofoil is clearly revealed in this 
figure. At a lift coefficient of 0.25 the L/D of the cambered hydro- 
foil is more than twice that of the flat plate. 

Also included in figure 31 are data taken from reference A on a 
ventilated modified flat plate of aspect ratio 1. This modified plate 
had an elliptical nose and a tapered trailing edge. The importance of 
providing a sharp leading edge on hydrofoils designed for use in cavity 

flow is shown by comparing the L/D of this modified flat plate with 

the L/D of the sharp-nosed flat plate of the present investigation. 

The rounded leading edge of the modified plate is subjected to a net 
positive pressure which is not balanced by similar pressures on the rear- 
ward portion of the plate. In addition, the lower surface of this plate 

presents an effective negative camber to the flow and thus does not develop 

as much lift as a truly flat surface. At small angles the drag of the 
rounded nose greatly influences the maximum L/D of the section. Several 
other investigations have noted the importance of a sharp leading edge on 
hydrofoils designed for operation in the supercavitating regime. (See 
refs. 18 and 19.) 

Center of pressure .- The center of pressure of the flat and cambered 
models calculated from the data of figures 22 and 23 are compared with 
theory in figure 32 . The theory from equation (64) is in good agreement 
with the experimental data for both models. Since the accuracy of the 
forces and moments on the flat plate is poor at small angles of attack 
(small total loads) the accuracy of the center of pressure from the data 
obtained on the flat plate at angles of incidence of 4° and 6° is doubtful. 


NONZERO CAVITATION NUMBER 


The nonzero-cavitation-number characteristics of the two models 
obtained at a 6 -inch depth of submersion are shown in figure 33 for angles 
of attack of 16°, 18°, and 20°. The solid data points are the lift coef- 
ficients obtained for vapor cavitation. The cavitation number corresponding 
to the condition tested was computed by using the water vapor pressure at 

the test temperature for the pressure within the cavity (that is 
p 0 - p„\ V 

o = — -I. For cavitation numbers less than about 0.7 the vapor pres- 


sure was the same as the measured cavity pressure. However, at cavitation 
numbers greater than 0 . 7 , the measured pressure was usually higher than 
the vapor pressure. Since it could not be determined with certainty 
whether the cavity pressure orifice was within the cavity, the vapor pres- 
sure was used to compute all cavitation numbers denoted by the solid data 
points. Also denoted in figure 33 are the approximate values of a at 
inception and the point at which the cavity length exceeded the chord. 
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These values are only estimates, since no effort was made to find the 
exact velocity at which these incidents occurred. Although the estimated 
values of are the same for both models, there may actually be some 

difference in the true points. 

The blank data points in figure 33 represent data obtained by intro- 
ducing air to the upper surface of the model and establishing a cavity. 

The cavitation number for this condition was computed by using the 
measured value of the pressure within the cavity. In this case the cavity 
pressure orifice was always well within the cavity formed. 

In figure 33 the curves are drawn through the computed value of 
for a = 0 as obtained from equation ( 50 ) . 

The agreement between the vapor- and air-cavity data confirms the 
use of the cavitation number as the significant parameter for correlating 
the characteristics of cavity flow. The similarity of the air and vapor 
cavities at nearly equal cavitation number is shown in figure 34. This 
dependence on the cavitation number is to be expected because the forces 
on the body are influenced only by the streamline curvatures and thus the 
pressure within the cavity; the type of gas present should have only a 
secondary influence. 

It was not possible to establish satisfactory air cavities at angles 
less than about 1 6 °. At low speeds where the air could reach the sepa- 
rated region near the leading edge the cavity upper surface was greatly 
disturbed by the force of the air jets. At higher velocities either the 
air could not reach far enough upstream to form a cavity or when a cavity 
was formed it did not cover the whole chord. If a greater quantity of 
air is supplied it is believed that satisfactory results can be obtained 
at angles less than 16 °. 

The use of a dynamic model of a high-speed aircraft will require 
simultaneous reproduction of both the full-scale cavitation and Froude 
numbers. Such an investigation is possible if an air-filled cavity with 
the proper cavitation number can be established on the model. 


CONCLUSIONS 


Conclusions based on the results obtained from the theoretical and 
experimental investigations of supercavitating flat and cambered hydro- 
foils may be summarized as follows: 

1. The theoretical expressions derived for the lift, drag, and 
center of pressure of supercavitating hydrofoils of arbitrary section 
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operating at zero cavitation number, and. finite aspect ratio and depth 
of submersion are in good enough agreement with the available experimental 
data to warrant their use in engineering calculations. 

2. The experimental ventilated force and moment coefficients of both 
models investigated at shallow depths of submersion were independent of 
speed. 

3 . Similar to the influence of camber on a fully wetted airfoil, the 
influence of camber in ventilated flow was to effectively increase the 
single of attack without appreciably changing the lift-curve slope. 

4. The lift-drag ratio of the cambered model near the design lift 
coefficient was more than twice that of the flat plate. 

5 . Comparison of the sharp-nosed flat plate with a rounded-nose flat 
plate showed the sharp-nosed section to be considerably superior. Thus, 
hydrofoils designed for operation in the supercavitating regime of flow 
should have sharp leading edges for best efficiency. 

6 . The cavitation number defines the flow similarity and lift of a 
supercavitating hydrofoil regardless of the type of gas (air or water 
vapor) within the cavity. Thus vapor cavities, normally obtained at high 
speeds can be simulated at lower speeds by establishing a cavity with 
forced air whose cavitation number is the same as that of the high-speed 
flow. 
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APPENDIX A 


THE TULIN -BURKART LINEAR THEORY FOR ZERO-CAVITATION-NUMBER 
FLOW APPLIED TO A CIRCULAR-ARC HYDROFOIL 


The problem is to find the effective increase in angle of attack due 
to the camber of a two-dimensional circular -arc hydrofoil operating at 
infinite depth and zero cavitation number by the method of reference 1, 
that is, to find the section coefficients Aq, Ap, and A 2 in the fol- 
lowing equation if a, measured from the chord line, is zero: 

C L = f (Ad' + A 1 " If) = f ( A o' + a c) 

By using the system of axes shown in the following sketch 



and the notation of reference 1 (barred symbols for the equivalent air- 
foil and unbarred for the hydrofoil) , the equation of the circular-arc 
hydrofoil is 


(x - c/2) 2 + (y + R cos y/2) 2 = R 2 (Al) 

il , . * - <=/ 2 - (A2) 

^ /r 2 - (x - c/2) 2 
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From reference 1 the slope of the equivalent airfoil is obtained from 
the equation 


*£(x) = Mx 2 ) 

dx dx' ' 


(A3) 


By noting that c = c , the slope of the equivalent airfoil section is, 
therefore, 


dy = 

dx 


x 2 - c g /2 


]r 2 - (x 2 - c 2 /2) 2 


(A4) 


If the conventional substitution of 


x = ~(l - cos 9) 


(A5) 


is made, equation (1*) after some manipulation may be written as 

.2 


but 


fly = 

dx 


i(l - cos 9) - 1 


1 

CVJ 1 

K] loj 

j|(l - cos 9) 2 - i| 

4 



(A6) 


R 2 _ R 2 
c^ c 2 V c /2 

4 4 


- IXf - 1- 


V S1„ | 


(A7) 


Therefore, 


dy 

dx 


1 

2 


cos 9-2 cos 9-1 


/ 1 -1 

f cos 2 9 - 2 cos 9-1) 

1 sin 2 Z 

l 2 ) 


(AS) 
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Since the linear theory is applicable only to small 7 and, 
therefore. 


1 >;> ( cos 2 9 - 2 cos 8 - l\ 

sin 2 1 \ 2 ) 

equation (A8) may be approximated as 

^ = - ^(cos^9 - 2 cos 9 - l) (A9) 


The necessary coefficients are then readily obtainable as 





A, = § f ^ cos 0 d9 = 2 

1 «J 0 to 2 

hr- = - f ^ cos 29 d9 = - 1 
xJq dx 8 


(A10) 


(All) 


(A12) 


An ~ 0 


(n > 2) (A13) 


The required effective increase in angle of attack due to small circular- 
arc camber is, therefore. 


a c = A 1 


A2 

2 


= 1.1 

2 2 


1 
8 > 



(A110 


It should be noted that for the reference line used in the analysis, 
A o' is not zero and positive lower surface pressures cannot possibly be 
realized near the leading edge unless the angle of attack is increased at 
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least to the point where a - 21 = 0 . Since = 0 (with n > 2) and. 

A]_ sin 9 + A 2 sin 29 is everywhere positive in the interval 0 ^ 9 = n, 

the condition a - 2. = 0 is sufficient to specify positive pressures 

a 

over the entire chord of the hydrofoil. 


A convenient way of treating the circular-arc section to make it 
comparable to other low-drag sections is to reorient its reference line 

so that a Q = that is, Aq' = 0. This orientation then corresponds 

to a = 0. At this design angle of attack = 0, a ° = 5) ^y < ^ ro_ 

dynamic efficiency Cj^Cp 

computed from linear theory is 21 , 

16 2C^ ^ 

low-drag hydrofoil selected in reference 1, ; — = -2- ]. 


of a two-dimensional circular-arc section as 
An 

— —X— which is almost as good as the 
/ Ct 


V C d 16 2C L; 


CONFIDENTIAL 


46 


CONFIDENTIAL 


NACA RM L57I16 


APPENDIX B 

SAMPLE CALCULATION OF LIFT, DRAG, AND CENTER OF PRESSURE 
OF A CAMBERED LIFTING SURFACE OPERATING AT FINITE 
DEPTH AND ZERO CAVITATION NUMBER 

A sample calculation of the lift, drag, and center of pressure is 
presented for a Tulin-Burkart section with A^ = 0.2 ^Cj^d. = 0*39 
having an aspect ratio of 1 and operating at an angle of attack of 12° 
and depth of submersion of 0.071 chord. 

LIFT COEFFICIENT 


Step 1 

For the Tulin-Burkart hydrofoil section 
tion (7), 



from equa- 


a r 


- A 1 - ¥ - TT ' °' 2 5 


Step 2 

From figure 11 at d/c = 0.071> 


a C,oo 


0.718 


Therefore , 

a c = (0.718) (0.25) = 0.18 
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Step 3 

By assuming that Cp q = 0.25 and r = 0.12 and by using equa- 
tion (L6) , 

cu = 1.12 ^ 2 - = 0.088 

x 7t 


and 


a + a c - aj_ 


12 

57.3 


+0.18 - 0.088 = O.3OI radian = 17.2° 


From figure 2 for d/c = 0.071 and a = 17.2°, it is found that m = I.63 
therefore, from reference 11, t = 0.12 as assumed. By using equa- 
tion (50) , 


C 


L, 1 " 


|(1.63) (0.301) 


cos 12° 
cos 17. 2° 


0.251 


This should check with the original assumption - if not, repeat step 3 
with better approximation for Cp and t. 


Step 1|- 

For the Tulin- Burkart section given by equation (12) with Aq = 0.2, 



Aq 

6 


Therefore, 


a 0 = tan _1 0.033 = I.92 0 


and 


a' = 12° + 1.92° = 13. 92° 
a c = 0.18 radian = 10. 3° 
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or 

a' + a c = 24. 22° 

By using equation (48) , 

C L ^ c = | 0.88 sin 2 (24. 22) cos 13.92° = 0.072 
Thus, the required lift coefficient is 

C L = C L x + C LfC = 0.251 + 0.072 = 0.323 


DRAG COEFFICIENT 


Step 1 

From step 3 in the lift-coefficient calculation. 


CL 4* CLq — CL^ — 17 • 2° 

For d/c = 0.071, m = I.63; and for d/c = <*>, m = 1.2. Therefore, from 
equation (6l) 


_ 1.63 * _ 

-TTl- 2 ' 13 


From figure 2 by using the a = 0 line with m = 2.13, it can be found 
that 

(d/c) e = 0.66 
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Step 2 

If a (the distance to airfoil center of pressure) is assumed to 
be 0.37* then from figure 9 for a = 0.37 and d/c = 0.66, the following 
values are found: 

B 0 = O.56 
B-l = 0.15 
b 2 = 0 

From equation (35) for Aq = a - = 0.121, A^_ = 0.2 and A 2 = -0.1. 

The values of Cq, Cj_, and C 2 may be determined as 

C 0 = 0.212 
C-l = 0.l8l 
C 2 = -0.1 

and from equation (29), a = O.367, which checks with the assumed value. 
If the resulting value of a differs enough from the assumed value to 
affect the values of the B coefficients, step 2 should be repeated. 


Step 3 

From equation (60) the value of e is determined as 

e = 0.0013 radian = 0.065° 

Step 4 

From equation ( 51) , 

C L, 1 = 0.251 C L ^ c = 0.072 C f = 0.004 

Therefore, 

Cp = 0.251 tan (12. 065°) + 0.072 tan (13.92°) + 0.004 = O.O75 
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CENTER OF PRESSURE 


Step 1 

From step 2 of the drag-coefficient calculation. 


C 0 = 0.212 Ci = 0.181 C 2 = -0.1 
Thus, from equation (62), 


x c.p.,l 


0.43c 


Step 2 

From equation (64) and hy using 


c L ,i = 0.251 


and 


c L,c = °-°7 2 


it is found that 


x c.p. = °- 4 5c 
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a 


£ 2 . = £2 £x for vapor cavitation 

q. q 




(b) Supercavitating or ventilated flow near the free surface. 
Figure 1.- Definition sketch. 
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Figure 2.- Green 's solution for lift-curve slope of a two-dimensional ventilated flat plate 
with approximate lines of constant ratio of depth of submersion to chord d/c. 
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Figure 3 • - Green's solution for effect of spray thickness on spray angle of two-dimensional 

flat plate. 
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Figure 4.- Two-dimensional theories for lift coefficient of 

cavitation number and infinite depth of 


circular -arc hydrofoil at zero 
submersion. 
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(a) Hydrofoil, Z-plane. 
Y 



(b) Equivalent airfoil, Z-plane. 

Figure 5 .- Linearized boundary conditions in hydrofoil and equivalent 

airfoil planes . 
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T 



Equivalent airfoil, S— .plane 


Figure 


6.- Linearized model for calculating effect 
sion on lift coefficient of flat plate (a 


of depth of submer- 

= 0 ). 



Figure 7.- Comparison of linearized solution with exact soultion for 
effect of depth of submersion on lift coefficient of a flat plate 
(a = 0) . 
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Equivalent airfoil } S- plane 


VJ1 
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Figure 8. - Linearized model for calculating effect of depth of submersion on vorticity distri- 
bution of equivalent airfoil section. 
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B 


0 



a, chords 


(a) B 0 . 



(*) B l . 

Figure 9.- The B coefficients. 
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(c) B 2 . 



(d-) B^. 



(e) B 4 . 

Figure 9 .- Concluded. 
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Figure 10.- Influence of depth of submersion on lift coefficient of various cambered sections 

operating at design angle of attack (a = 0) . 
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Figure 11.- Influence of depth of submersion on effective angle of attack due to camber <x c of 
various cambered sections operating at design angle of attack (a = 0) . 
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(a) 5/c = 00 . 


ii 

c 


Free surface 


v 


(b) 6/c = 0. 



Figure 12.- Location and direction of stagnation line on flat plate 

normal to direction of motion. 
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(b) Airfoil, Z = -\J~Z plane. 

Figure 13 .- Hydrofoil and equivalent airfoil at large angles of attack 

and small depths of submersion. 
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(a) Cambered model. 

Figure l4. - Model configurations. 
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(b) Flat-plate model. 
Figure 14.- Concluded. 
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L- 57- 27^5 

Figure 15 .- Test set-up showing cambered hydrofoil with aspect ratio 
of 1 and balance attached to towing carriage. 
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Lift 



->-• Direction 

of motion 


Figure 16.- System of axes. 




Figure 17.- Schematic drawing of spray-thickness gauge. 
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V = 20 fps 



V = 60 fps 


(a) a 


V = 40 fps 






V = 20 fps 




V = 30 fps 


(b) a = 10 °. 


V = 40 fps 


L- 57 - 2740 


Figure 18.- Flow about the flat lifting surface, depth of submersion is 

0.5 inch. 
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V = 30 fps 


V = 40 fps 


(c) a = 8°. Air supplied at 0.012 pound per second for V = 10, 20, 

30 feet per second. 


Figure 18.- Concluded. 



V = 20 fps 


V = 60 fps 


V = 40 fps 



V = 80 fps 



(a) a = 10°. 


L- 57-2741 


Figure 19.- Flow about cambered lifting surface, depth of submersion is 

0.5 inch. 
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V = 20 fps 



V = 40 fps 



V = 80 fps 


Cb) a = 





V = 40 fps 

16°. L- 57- 2742 


Figure 19 . - Concluded. 
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Figure 20.- Vortex ventilation speed for 0.5 inch depth of submersion. 
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(a) Flat plates. 



Figure 21.- Theoretical location of cavity upper surface for two- 
dimensional flat and cambered lifting surfaces, (infinite depth; 
a = 0 .) 
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(a) Lift. 


Figure 22.- Characteristics of ventilated flat lifting surface. 
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(b) Drag. 


Figure 22.- Continued. 
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(c) 


Pitching moment about leading edge. 


Figure 22.- Concluded. 
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(a) Lift. 

Figure 23.- Characteristics of ventilated cambered lifting surface. 


CONFIDENTIAL 


NACA RM L57I16 


CONFIDENTIAL 


79 


440 


400 


360 


Depth, 

a, 

in. 

deg 

O 1/2 

20 

□ 1/2 

IS 

0 1/2 

16 

A 1/2 

14 

k. 1/2 

12 

& 1/2 

10 

O 1/2 

8 

a 0 

10 

a 0 

8 

O' 0 

6 


Ventilation 

Vortex 
Vortex 
Vortex 
Vortex 
Vortex 

Vortex, wire on L.E. 
Forced air, wire on L.E. 
Vortex 
Vortex 



40 50 

Speed, fps 


(b) Drag. 

Figure 2J.- Continued. 
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(c) Pitching moment about leading edge. 


Figure 23.- Concluded. 
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(a) Lift. 


Figure 2k. - Characteristics of ventilated flat plate at 
for depths of submersion of 0, 1.0, 1.5, and 2.0 


a = 16 ° and 20° 
inches . 
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(b) Drag. 

Figure 24.- Continued. 
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(c) Pitching moment about leading edge. 
Figure 24.- Concluded. 
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Figure 25.- Effect of depth of submersion on flat lifting surface spray thickness with aspect 

ratio of 1. 
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Angle of spray with the horizontal. 



Figure 26.- Effect of depth of submersion on flat lifting surface spray angle with aspect ratio 

of 1. 
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C □ 1/2 Experiment, vortex ventilation 
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Figure 27.- Comparison of theoretical and experimental ventilated lift coefficients for flat 

and cambered lifting surfaces with aspect ratio of 1. 
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Figure 28.- Comparison of theoretical and experimental lift coefficients for flat lifting sur- 
face as affected by depth of submersion, a = l6° and 20°; aspect ratio of 1. 
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O □ <0> A Experiment 
Theory 



(a) Cambered section, a = 9°- 



(b) Flat plate. 

Figure 29.- Comparison of present theory with experimental data of 
Fuller obtained on hydrofoils having an aspect ratio of 2. 
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Figure 50.- Comparison of theoretical and experimental ventilated drag coefficients for flat 

and cambered lifting surfaces with aspect ratio of 1. 
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Figure 31 . - Comparison of theoretical and experimental ventilated lift- 
drag ratios for flat and cambered lifting surfaces with aspect ratio 
of 1. 
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Figure 32.- Location of center of pressure of flat and cambered lifting surfaces with aspect H 

ratio of 1. 
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(a) a = l6°. 



(b) a = 18°. 



(c) a = 20°. 


Figure 35.- Effect of cavitation number on lift coefficients of flat 
and cambered lifting surfaces. Depth of submersion is 6 inches; 
aspect ratio of 1. 
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( a) Flat plate . 


L- 57-27^3 


Figure jk.- Comparison of air- and vapor-filled cavities for depth of 
submersion of 6 inches and a of 20°. 
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Air, 0 = 0.48, V = 20 fps 



Vapor, cr = 0.35> V = 80 fps 

(b) Cambered section. L-57-27^ 

Figure 3b. - Concluded. 
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